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Abstract — In a recent work, Nazer and Gastpar proposed 
the Compute-and-Forward strategy as a physical-layer network 
coding scheme. They described a code structure based on nested 
lattices whose algebraic structure makes the scheme reliable and 
efficient. In this work, we consider the implementation of their 
scheme for real Gaussian channels and one dimensional lattices. 
We relate the maximization of the transmission rate to the lattice 
shortest vector problem. We explicit, in this case, the maximum 
likelihood criterion and show that it can be implemented by using 
an Inhomogeneous Diophantine Approximation algorithm. 



I. Introduction 

In [1], Zhang et al. introduced the physical-layer network 
coding concept (PNC) in order to turn the broadcast property 
of the wireless channel into a capacity boosting advantage. In- 
stead of considering the interference as a nuisance, each relay 
converts an interfering signal into a combination of simulta- 
neously transmitted codewords. PNC concept has received a 
particular interest in the last years because it provides means 
of embracing interference and improving network capacity. 

In a recent work [2], Nazer and Gastpar proposed a new 
physical-layer network coding scheme. Their strategy, called 
compute-and-forward (CF), exploits interference to obtain 
higher end-to-end transmission rates between users in a net- 
work. The relays are required to decode noiseless linear 
equations of the transmitted messages using the noisy linear 
combination provided by the channel. The destination, given 
enough linear combinations, can solve the linear system for 
its desired messages. This strategy is based on the use of 
structured codes, particularly nested lattice codes to ensure 
that integer combinations of codewords are themselves code- 
words. The authors demonstrated its asymptotic gain using 
information-theoretic tools. 

The authors in [3] followed the framework of Nazer and 
Gastpar and showed the potential of the compute-and-forward 
protocol using an algebraic approach. They related the Nazer- 
Gastaper's approach to the theorem of finitely generated 
modules over a principle ideal domain (PID). They gave 
sufficient condition for lattice partitions to have a vector 
space structure which is a desirable property to make them 
well suited for physical-layer network coding. Then, they 
generalized the code construction and developed encoding and 
decoding methods. 

In [4], the authors proved that the lattice implementation 
of compute-and-forward as proposed by Nazer and Gastpar 
suffers from a loss in number of achieved degrees of freedom. 
They proposed a different implementation consisting of a 
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Fig. 1. System model: 2 sources and one relay. 

modulation scheme and an outer code and showed that it 
achieves full degrees of freedom as if full cooperation among 
transmitters and among relays was permitted. In their scheme, 
the channel coefficients are known throughout the network. In 
[5], the authors designed a modulation/coding scheme inspired 
by the compute-and-forward protocol for the wireless two-way 
relaying channel. 

In this work, we consider the practical aspects of the 
compute-and-forward protocol. We implement the protocol 
described by Nazer and Gastpar. We explain how to obtain the 
integer coefficients that maximize the rate. We also propose 
a decoding technique based on maximum likelihood. Finally, 
we show some simulation results. All the practical aspects are 
demonstrated here for one-dimensional real constellations. 

II. System Model and Assumptions 

In our model, we consider one relay receiving messages 
from two sources S\ and S 2 and transmitting a linear combi- 
nation of these two messages, as described in Figure Q] The 
relay observes a noisy linear combination of the transmitted 
signals through the channel. Received signal at the relay is 
expressed as, 

y = hixi + h 2 x 2 + Z. (1) 

The relay searches for the integer coefficient vector a = 
[ai a 2 ] T that maximizes the transmission rate. It then decodes 
a noiseless linear combination of the transmitted signals, 



xr = a\X\ + a 2 x 2 , 



(2) 



and retransmits it to the destination or another relay. We con- 
sider a real- valued channel model with real inputs and outputs. 
The channel coefficients hi and h 2 are real, i.i.d. Gaussian, 
hi ~ AT(0, 1). z is Gaussian, zero mean, with variance a 2 = 1 
(z ~ A/"(0, 1)). Let h = [hi h 2 ] T denotes the vector of 
channel coefficients. Source symbols xi are integers and verify 
\xi\ < s m , i.e., Xi £ S = {-s m , -s m + 1, ■ • ■ , s m }. Si and 
S 2 transmit xi and x 2 , respectively. Both sources have no 
channel side information (CSI). CSI is only available at the 
relay. 



III. COMPUTE-AND-FORWARD 

In what follows, we use the expression of the computation 
rate R comp given by Nazer and Gastpar [2] in order to find 
a vector a maximizing it. We show that the maximization 
of i? C omp is equivalent to the search of a shortest vector in 
a lattice. Then, based on the likelihood expression, we show 
that decoding is equivalent to an Inhomogeneous Diophantine 
Approximation. 

A. Achievable Computation Rate 

The primary goal of the decode-and-forward is to enable 
higher achievable rates across the network. Nazer and Gastpar 
showed that the relays can recover any set of linear equations 
with coefficient vector a as long as the message rates are less 
than the computation rate 



Rcom P {h,a) =log 
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(3) 



where this rate is achievable by scaling the received signal 
by the MMSE coefficient [2]. We are interested in finding the 
coefficient vector with the highest computation rate. This is 
given in the following theorem. The result is obtained for a 
relay combining N symbols and for complex- valued channels. 



Theorem 1: For a given h s C (resp. R ), R comp (h,a) 
is maximized by choosing a € %[i] N (resp. Z ) as 



a 



arffmin (a^Ga) 



where 
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(4) 



(5) 



H = [Hij], H %3 = hih*, 1 < i, j < N and f is for the 
Hermitian transpose (resp. the regular transpose). 

Proof: Maximizing R comp (h,a) is equivalent to the 
following minimization 
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We can write 



As H 



SNR||ft|h|ar-SNR 



\h'a\ 2 = 2_[ hih*a*aj 



h)a\ 2 } 



(6) 



[Hij], H i3 = hih*, 1 < i,j < N, it follows that 
J2hih*a*a,j = a'' Ha. Using these notations, we can write 

© as 

SNR 



(l + SNRIIhlHminat 



tH 



(8) 



SNR 



l + SNR||h|| 2 

H has N strictly positive eigenvalues. It 

is then positive definite. Now, the problem is reduced to the 
minimization of a^Ga. ■ 

Proposition 1: Searching for the vector a that minimizes 
Equation © of theorem Q] is equivalent to a "Shortest Vector" 
problem for the lattice A whose Gram matrix is G. 

Proof: As G is a definite positive hermitian (resp. 
symmetric) matrix, it is the Gram matrix of a lattice A. This 



lattice is either a Z[i]— lattice in the complex case, or a Z— 
lattice in the real case. Then, the minimization problem in 
theorem Q] is equivalent to find a non zero vector in A with 
shortest length. ■ 

Algorithms for solving this problem are given in [6], The 
best known one is the Fincke-Pohst algorithm [7]. 

B. Recovering Linear Equations 

The relay aims to decode a linear equation of the transmitted 
messages and passes it to the destination or another relay. After 
calculating the vector a as in ©, the relay recovers a linear 
combination of the transmitted signal x\ and x%. We rewrite 
the received signal at the relay in the following form 



y = X + £,lXx +&X2 + z 



(9) 



where A is an integer, £j = hi — cii and z is the additive 
white noise. The recovered linear equation A = a\X\ + 0,2X1 
is a linear Diophantine equation. This equation admits the 
following solutions. 

C. Solution of the Linear Diophantine Equation 

If A is a multiple of the greatest common divisor (gcd) of a\ 
and <i2, then the Diophantine equation has an infinite number 
of solutions. The Extended Euclid Algorithm allows to exhibit 
a particular solution (u\,U2) to a\X\ +022:2 = 9 [9]. The set 
of all solutions is obtained as follows 




9 9 



(10) 



9 9 

a\ A 0,2 is the gcd of a\ and 0,2, k E Z. 



D. Decoding Metric 

The Maximum Likelihood decoder maximizes p(y/X) over 
all possible values of A. The conditional probability p(y/X) 
can be expressed as, 



p(y/x) 



^2 p(v/ x i> x 2)p(xi,x 2 ) (ii) 



(xi,x 2 ) 



(') where 



p(y/xi,x 2 ) oc exp 



{y - hixi - h 2 x 2 )' 2 



2a 2 



(12) 



and x\, X2 are (a priori) equiprobable and given by dTOb . The 
decoding rule is now to find, 



+00 



A = argmaxg(A) := >J exp 



fc=— 00 



(y-/3X + kaf 
2^ 



(13) 



where j3 — | (hiUi + h 2 u 2 ), a — A (h 2 ai - h 1 a 2 ). 

In [8], it has been proved that, for A 6 R, q(X) achieves its 
maximum for 

i.e. for all values of A such that y — f3X + ka = 0. Since 
we want to maximize g(X) for A 6 Z, the solution is given 
by the integer- valued couple (A, k) minimizing \y — (3X + ka\. 



Thus, since x\ , X2 G S which is a finite subset of Z and verify 
Equation ( [Tol l, we state a new minimization problem which is 
equivalent to (Qjji, 
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(14) 



The problem is therefore equivalent to the minimization of 

F(k,\) = \ka'-\ + y'\ (15) 

a' = a//3 and y' = y//3. The minimization is called Inhomo- 
geneous Diophantine Approximation in the absolute sense. It 
consists of finding the best approximation of a real number a' 
by a rational number X/k, k 6 N, given an additional real shift 
y 1 ', while keeping the denominator k as small as possible. In 
the general settings for such problems, an error approximation 
function F(k, A) is set and it is stated that a rational number 
X/k is the Best Diophantine Approximation if, for all other 
rational numbers A'/fc' 



k' < k=>F{k',X') >F(k,X). 



(16) 



In our case, in addition to the error approximation function, 
limits are imposed by the finite constellation S to which 
the transmitted symbols belong. The algorithms used to find 
the best Diophantine approximations of real numbers are in 
general simple and easy to implement. The best known one is 
the Cassel's algorithm [10]. In [11], the authors develop and 
compare several ones. 

IV. Numerical Results 

In the simulations, the set of symbols is of the form S = 
{— s m , . . . , s m }. We consider two sources transmitting x\ and 
X2, and one relay recovering a linear equation of x\ and X2 
with integer coefficients. 

At first, based on its CSI, the relay finds the vector a as 
the shortest vector described in theorem Q] Then, the relay 
finds a particular solution of the linear Diophantine equation 
aiXi+a,2X2 — g using the Extended Euclid algorithm. Finally, 
the relay searches for the couple (k, A) which gives the best 
inhomogeneous Diophantine approximation by minimizing the 
function F defined in ( [TBI ). 




Fig. 2. p(y/X) for h = [-1.274 0.602] T , a = [2 - 1] T , SNR = 40<£B, 
x\ = —2 and X2 = 3. p(j//A) is maximized for one value, A = —7 in the 
left subfigure while it is maximized for several values of A in the right one. 

In Figure [3] we show the error probability of our system 
for three different constellations S, defined by s m = 5, 7, 10, 
respectively. For s m = 5 or less, the diversity order of the 
system is 1 for real entries (which would correspond to a 
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Fig. 3. Error Probability using the Inhomogeneous Diophantine Approxima- 
tion versus decoding both symbols x\ and X2- 



diversity order equal to 2 with complex symbols). For s m > 6, 
the diversity order collapses to 1/2. This is due to the fact that 
p(y/X) is constant, as a function of A, on a bigger interval 
giving rise to ambiguities as shown in Figure [2] Still in Figure 
|3] we plotted the error probability for when the relay decodes 
both symbols x\ and X2- The diversity order in this case is 1/2 
for all values of s m . For the case of complex- valued channels 
and symbols, we expect a doubled value of all the diversity 
orders. 

V. Conclusion 

In this paper, we considered the Compute-and-Forward 
scheme with real-valued channels. We provided a method for 
maximizing the transmission rate and developed a decoding 
strategy. Numerical results showed the performance of our 
decoding method. We believe that it is a first step towards 
a rich and fruitful multidimensional approach. 
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